Continuum in the Excitation Spectrum of the S = 1 Compound CsNiCl3 
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Recent neutron scattering experiments on CsNiCl3 reveal some features which are not well de- 
scribed by the nonlinear a model nor by numerical simulations on isolated S=l spin chains. In 
particular, in real systems the intensity of the continuum of multiparticle excitations, at T = 6 K, 
is about 5 times greater than predicted. Also the gap is slightly higher and the correlation length 
is smaller. We propose a theoretical scenario where the interchain interaction is approximated by a 
staggered magnetic field, yielding to a correct prediction of the observed quantities. 

PACS numbers: 75.10. Jm, 75.45.+j, 75.50.-y 



About twenty years after the famous "Haldane 
conjecture" 0], quantum spin chains with Heisenberg an- 
tiferromagnetic interactions are still attracting much ex- 
perimental and theoretical interest. Half-integer spin 
chains have no spin gap, and hence are critical at T = 
with algebraic decay of correlations. Integer spin chains 
have a gap in the excitation spectrum and the ground 
state is characterized by a finite correlation length. 

For S=l chains, the excitation spectrum is dominated 
by a well-defined magnon triplet carrying spin 1. These 
excitations have been observed in several neutron scatter- 
ing experiments and in the case of NENP, they constitute 
a large part of the total spectrum. 

However, recent experimental measurements on 
CsNiCl3 reveal a significant spectral weight in the in- 
coherent multiparticle continuum above the coherent 
magnon peak. At 6 K the integrated intensity of the 
continuum around the antiferromagnetic point is about 
9(2)% of the total spectral weight^). This result is 
considerably larger than the 1-3% weight predicted by 
numerical simulations^ and the 0(3) non-linear sigma 
model for the 3-particle continuum The effects of 
the coupling between chains has been considered in the 
framework of RPA, but no significative increase of the 
continuum has been found Q. On the other hand, Ma- 
jorana fermion theory yields a prominent three-particle- 
scattering continuum but this would require strong 
biquadratic exchange interactions of the form (Sj • S^+i) 2 
which are not in accordance with the measurements on 
CsNiCl 3 §. 

Furthermore, some measured quantities at T = 6 K 
show appreciable deviations from analytical and nu- 
merical predictions at finite temperature on isolated 
chains H ||, i.e. the measured gap A is higher and the 
correlation length £ is shorter than the corresponding 
theoretical values, breaking also the single mode approx- 
imation (SMA) correspondence A£ = c, where c is the 
spin-wave velocity. 



The model Hamiltonian proposed for CsNiCLj is 
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where the coupling J = 2.28 meV along the c-axis is 
much stronger than the exchange interaction in the basal 
plane J' = 0.044 meV@. The single-ion anisotropy is es- 
timated to be D k 4 fieV, small enough that CsNiCi3 
may be considered purely isotropic, hence we put D = 
from now on. Below Tn = 4.85 K the system under- 
goes 3D long-range ordering, caused by the interchain 
interaction. For T > Tn the system is regarded as one- 
dimensional and commonly considered to be a good real- 
ization of a single isotropic spin chain, although there is 
a sizable dispersion perpendicular to the chain direction 
due to the interchain coupling. 

Here comes our main assumption: 
In some range of intermediate temperatures above the 
ordering one, T/v < T < T* , the interchain coupling 
induces a (short-range) three-dimensional antiferromag- 
netic order which varies on a length scale which is 
greater than the correlation length £id of the single chain. 

In other words, a little above Tjv, even if the total 
staggered magnetization is zero, three-dimensional ef- 
fects are still important and the system is arranged in 
large domains with non-zero magnetization. Thus, in a 
mean-field picture, a single chain experiences a staggered 
magnetic field generated by the neighboring chains. No 
matter if the field is not constant along the whole chain, 
because it varies so slowly that it may be considered con- 
stant, since the ID chains are gapped and short-ranged. 

Consistently, we suppose the existence of a tempera- 
ture T* above which the antiferromagnetic order is so 
reduced that the ID system may be viewed as isolated. 
To support this hypothesis we mention that in the real 
compound CsNiCi3, some experimentalists Q found that 
for T > 12 K the measurements turn out to be in accor- 
dance with the theoretical predictions on isolated chains. 

In the mean-field approximation the interchain inter- 
action is treated as a staggered magnetic field that takes 
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into account the effect of the neighboring chains on the 
single ID system 



plane 

(ij) 



H s j2(-ys?, (2) 



where the staggered field (we choose it along the z-axis) 
is evaluated self-consistently by means of the fixed-point 
equation B 



H s = z c J' m s (H s 



(3) 



with z c the coordination number in the basal plane and 
m s (H s ) is the magnetization curve of a single chain in 
presence of the external staggered field H s . In our cal- 
culation we fix the energy scale J = 1, the interchain 
coupling J 1 = 0.02 as close as possible to the experimen- 
tal value for CSMCI3 and we put z c = 3 because the 
lattice in the basal plane is triangular. 

So, the Hamiltonian we are going to study is 



(4) 



Following the Haldane mapping we represent the 
local spin operators as Si » S(— l) l nj + lj, with nf = 1, 
where represents the slowly-varying local staggered 
magnetization and lj is the local generator of angular mo- 
mentum. In this framework, the Zeeman term of Eq. (Jij) 
becomes essentially a linear shift in the n-field along the 
^-direction, (-l) 2 iJ s Sf « S H s n\. 

Going then to the continuum limit and integrating out 
the fluctuation field 1 we obtain the 0(3) NLcM plus a 
linear term]To| 
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where g = 2/S and c = 2JS. A Lagrange multiplier 
A(a;,r) has been introduced to implement the local con- 
straint n 2 (x,r) = 1. 

The experiments were performed at temperatures 
small enough if compared with the exchange interac- 
tion J along the chain, so we mainly consider the zero- 
temperature case arguing that a calculation at finite 
temperature j{| will give simply a small correction in the 
calculated quantities. 

For a constant staggered field, the associated saddle 
point will correspond to a constant value of A and the 
self-consistency equation at T = is 



^ln{Ae+ v / l + (A0 2 } 



= 1-1^ 
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where £ is the correlation length and A is the momentum 
cutoff of the theory. In the passage to the continuum 



limit we have lost the correct renormalization of the pa- 
rameters. In order to restore them, we fix the zero-field 
values of A = 0.41048 J, £ = 6.03 (which correspond 
to c = £A = 2.48 J), known from the DMRG studies. 
The NLctM coupling has been left to its analytical value 
g = 2. As the field is varied, the cutoff A is held fixed 
to its zero-field value A = 0.2072 and we use Eq.^) to 
determine the value of £ = £,(H S ). 

The magnetization per site is given by 



m s {H s ) = S(n z ) = 9 —e{H s 
c 



(7) 



Putting it back in Eq.(||), we find a stable solution for 
H s = 0.027, corresponding to m s = 0.45, £ = 4.54 and 
Ay = c£ _1 = 0.545. It is remarkable that even for a 
small staggered field the response of the system is strong, 
generating an appreciable magnetization, meanwhile the 
gap renormalizes upward and the correlation length is 
getting shorter. 

The staggered field breaks explicitly the symmetry and 
the quasi-particle propagator splits in one longitudinal 
channel (i.e. in the direction of the field) and two trans- 
verse channels. The transverse propagator is given by 
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and has simple poles at to — ±s{q) with 



e(q) = Jc 2 q 2 + A 2 T , 



(9) 



where q = corresponds to the antiferromagnetic point. 
In the transverse channel the theory is purely bosonic 
and the spectral weight is fully exhausted by this magnon 
excitation. In this case, the dynamical structure factor 
Sriq,^) = Im Gt (q, <jj)/it is simply given by 
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{S(uj-e(q))-S(u; + e(q))} (10) 



On the other hand, the longitudinal propagator must 
be evaluated more carefully considering its connected 
part and the calculation m/m gives 



G L (q,w) = G T {q,uj) 
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where T(q, tv) is the Fourier transform of the product of 



two propagators 

r(x-x') = ^jGr(x-x')Gr(x'- X ) . (12) 

From the study of the analytic structure of the longitudi- 
nal propagator it is possible to calculate the longitudinal 
gap solving the following equation for A^ [jy] 
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where Ti (0, A/,) is the real part of the "polarization bub- 
ble" (in q = and u) = A/,) and may be written as 
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For the field strength that we have estimated for CsNiCbj, 
we find the value Al = 0.779, which is greater than 
Ay, as we expected. It is important to stress that the 
presence of the continuum makes inapplicable the SMA, 
which establishes the relation \l = Sgc/ A| between the 
longitudinal susceptibility xl — dm s /dH s and the gap. 

For u> < 2At the dynamical structure factor in the 
longitudinal channel has well-defined poles corresponding 
to single particle excitations 



S L (q,u) = 7 



{Siuj-SL^-Siuj + eUq))} 



(15) 



where the prefactor 7 is less than unity and gives 
the reduction of the quasi-particle weight. From the 
calculation |p"i"f , it turns out that 7 is a decreasing func- 
tion of the staggered field, as displayed in Fig.|l|. 
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FIG. 1: The relative quasi-particle weight of the pole of the 
longitudinal propagator as a function of the staggered field, 
at q = 0. 
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FIG. 2: Continuum contribution to the spin dynamical struc- 
ture factor at the antiferromagnetic point, for H a = 0.027 
(solid line), H s = 0.020 (dotted line) and H s = 0.010 (dashed 
line). The integrated intensity grows as the staggered field 
increases. 



As the field increases, the spectral weight that is lost 
from the pole gets transferred to the multi-particle con- 
tinuum, as required by the first moment sum rulepd|] 
of the dynamical structure factor. We can observe this 
phenomenon by calculating the continuum contribution 
to the dynamical structure factor for some values of the 
staggered field. In Fig.|| the quantity S{q,io) is plotted 
as a function of uj at the antiferromagnetic point, for 
H a = 0.027, H s = 0.020 and H s = 0.010 (dashed line). 
The integrated intensity grows as the staggered field in- 
creases. We also notice that this contribution starts at 
the two-magnon threshold ui > 2At(H s ), consistently 
with experimental data|^, |l^| in which the continuum 
starts well before the three magnon contribution as pre- 
dicted by the NLctM for a single chain 

At H s = 0.027 the value of the prefactor is 7 = 0.746, 
hence using the first moment sum rule we find that 
the relative weight of the incoherent contribution to the 
structure factor is 
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.5% 



(16) 



which is a considerable amount, very close to the exper- 
imental measure of 9(2)%. 

The present theory predicts the existence of a transver- 
sal mode and a longitudinal one with different energies. 
This fact would be consistent with the presence of two 
well-defined peaks in the inelastic scattering measure- 
ments with unpolarized neutrons. However, in the ex- 
periments only one peak is observed. Remembering that 
a similar problem arise in the case of i?2BaM05^2|, we 
argue that the longitudinal mode could be masked by its 
proximity to the transversal one, by its reduced weight 
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(less than 1 /2 of the transverse weight) and by the wide 
broadening. 

In conclusion, we have proposed a theoretical mech- 
anism to explicate the anomalous multiparticle excita- 
tion weight observed in recent neutron scattering experi- 
ments. This quantity turns out to be much greater than 
the theoretical prediction on isolated 5=1 chains. We 
assume the existence of a temperature region above Tjv, 
where the 3D effects are not at all negligible, although 
the long range order is destroyed, and we take into ac- 
count of the interaction with the neighboring chains by 
means of a staggered magnetic field. We found that, 
even in the framework of a simple mean-field theory, this 
scenario may account for the generation of a sizable con- 
tinuum weight, as well as an upward gap renormalization 
and a reduction of the correlation length, breaking also 
the SMA correspondence between these two quantities. 
Further and more detailed analytical calculation and nu- 
merical tests will be treated in a future paper. 
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